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1 Introduction 

A system which is part of a larger system and interacts with its other parts cannot be in 
a pure quantum state and must be described by a density matrix [19]. The calculation 
of the density matrix involves taking the trace over all those degrees of freedom of the 
larger system which do not belong to the subsystem we are interested in and, in the 
thermodynamic limit, when the subsystem stays finite, but the large system becomes 
infinitely large, is usually rather hard. One of the few examples where a density matrix 
could be calculated for an interacting system is the antiferromagnetic spin-^ XXZ 
chain with Hamiltonian 

Hxxz = jt i^U^j + o]_,oy+A{a)_,<fj - 1)) (1) 

acting on the tensor product of spaces of states of L spins ^. The o", a = in (1) 
are the Pauli matrices and / > and A > —1 are two real parameters, the exchange 
interaction and the exchange anisotropy. 

In [8,9, 14] integral formulae for the zero temperature density matrix elements of 
a segment of length m of the infinite chain were obtained. Since the expectation value 
of any operator acting on a segment of length m can be expressed in terms of the 
density matrix elements (see equation (13) below), the density matrix enables one, in 
particular, to calculate the correlations of local observables (for recent developments 
see [13]). This is the reason why the density matrix elements were called 'elementary 
blocks of correlation functions' in [8,9, 14] and also is the actual reason why we got 
interested in the subject. In fact, much recent progress in the calculation of short 
range correlations for the XXZ chain [2, 3, 5, 10, 1 1, 20, 26] originates in the integral 
representation of the density matrix obtained in [8, 9, 14]. 

Below we generalize the formulae first obtained in [8,9, 14] to finite temperatures. 
It turns out that the same means that where successfully applied in the calculation of 
a generating function of the S^-S^ correlation functions at finite temperatures in [7] 
also work for the density matrix. Namely we may combine algebraic Bethe ansatz 
techniques for the calculation of matrix elements [14, 17, 18,21] with the quantum 
transfer matrix approach to the thermodynamics of quantum spin chains [15, 16,22,25]. 



2 The quantum transfer matrix 



All properties of the XXZ chain can be derived from the well-known trigonometric 
solution 

/I 0\ 

b{l,{j) c{X,/j) 
c(\//) b{'k,ij) 
\0 1/ 



(2) 
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of the Yang-Baxter equation, where 

sh(X-A/) s. sh(ri) 

H^^M) = -rpi : — r , c(k,IV = -rp{ : — r • (3) 

sh(X-A' + ri) sh(X-^ + ri) 

This -matrix not only generates the Hamiltonian (1), 

L 



Hxxz = 2ysh(Ti) d^{PR)j_ij{X,0) 

7=1 



(4) 



(P is the permutation matrix, ch(ri) = A in (1)), but also a related auxiliary vertex 
model whose partition function in a certain 'Trotter limit' is equal to the partition 
function of the XXZ Hamiltonian. Details of this construction in a notation also suit- 
able for our present purpose were reviewed in our paper [7]. Here we shall repeat 
only the most important formulae as far as they are needed to understand the notation 
below. We define the monodromy matrix of the auxiliary vertex model as 



\C{X) D(X)J. 

where ti means transposition with respect to the first space, and the parameter (3 is 
inversely proportional to the temperature T, 

3=^. (6) 

The monodromy matrix (5) acts in the tensor product of an auxiliary space j and 
'quantum spaces' 1, . . . ,A^. It generates a representation of the Yang-Baxter algebra 
with /?-matrix (2). The corresponding transfer matrix 

t{X)=trjTj{X), (7) 

called the quantum transfer matrix, defines a vertex model, whose partition function in 
the Trotter limit AT — > oo is equal to the partition function Zl = tre.xp{—Hxxz/T) of the 
XXZ chain of length L, 

oo 

Z, = ^im tri...^(r2^^(0))^ = £ A^(0) . (8) 

Here A„(X) denotes the nth eigenvalue of the quantum transfer matrix. Note that there 
is a unique real eigenvalue Ao(X) with largest modulus which dominates the partition 
function in the thermodynamic limit, when L goes to infinity. This single leading 
eigenvalue determines the bulk thermodynamics of the XXZ chain. We showed in [7] 
that the corresponding eigenvector determines the state of thermodynamic equilibrium 
completely. It fixes all finite temperature correlation functions (compare [23, 24]). 
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Due to the conservation of the z-component of the total spin 

thermodynamics and finite temperature correlation functions can still be treated within 
the quantum transfer matrix approach if the system is exposed to an external magnetic 
field h in z-direction [7]. The external field is properly taken into account by applying 
a twist to the monodromy matrix (5), 

r(X)^r(X)(^7^_o,,). (10) 



3 Integral formula for density matrix elements 

We shall use the notation 

for the gl(2) standard basis. The canonical embedding of these matrices into the space 
of operators on the space of states (C^)®^ of the Hamiltonian (1) will be denoted by 
e^p , a, P = 1 , 2, J = 1 , . . . , L. Then every operator A\^_,_^rn that acts on sites 1 to m of 
the spin chain can be expanded as 

Al,...,m =^ai...oc1^1pJ • • -^mp™ ' (12) 

where implicit summation over the Greek indices is implied. The thermal average of 
such type of operator is 

{Ai,...,m)T,h=At-k{eil\ ■ ■■em;:)r,h ■ (13) 

Thus, it is sufficient to calculate the expectation values . ■^^m^'")J•f^ which define 
the matrix elements of the density matrix of a chain segment of length m. 

Following [7] we can calculate the general density matrix element as a limit of an 
appropriately defined inhomogeneous finite Trotter number approximant, 

where 

for Uj.^k =1,2, and = . . .B(k^^2)\^) is the eigenstate corresponding to 
the leading eigenvalue Ao(A,) which is parameterized by a specific set of Bethe roots 
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j 


1 2 3 4 5 6 7 8 


Pi 


i i T T i T i i 

i T i i i T T i 


a; 
P7 


3 4 6 

1 3 4 5 8 


a; 
P7 


6 4 3 

1 3 4 5 8 




1 3 3 4 4 5 6 8 


yr 


3 5 7 




1 2 4 6 8 



Table 1: Example for the definition of the sequences {clJ), (P^- ), (oc^) and (P^ ). 
The pattern corresponds to the string D(^i)C(^2)5(^3)5(^4)^(^5M(^6)C(^7)^(^8) 
of monodromy matrix elements, m = 8, |a+| =3, |(3 | =5. The sequences {jj) and 
(jf) are defined in section 4 below where they are needed. 

{X} = The complex inhomogeneity parameters 'regularize' the expres- 

sion in the numerator on the right hand side of (15). Moreover, at least in the zero 
temperature case, it turned out to be useful and interesting to study the dependence of 
Dat on these parameters [1,4]. 

We may identify values 1 of the indices with the symbol |, representing an up- 
spin, and values 2 with |, representing a down spin. Then the upper and lower indices 
(o'n)™=i and (|3n)™=i in e.g. (15) may be visualized as sequences of up- and down- 
spins. We shall denote the position n of the jth up-spin in the sequence (aw)^=i by 
a+ and the number of up-spins in (a„)^^j by |a+|. Similarly we define Pj as the 
position of the jth down-spin in (Pn)J^=i and denote the number of down-spins in 

(Pn)n=i by I PI- This yields two sequences (at)l^''l and {^J)%1 which we reorder 
and shift by the prescription at = aj^+|_^._^j, 7 = 1, . . . , |a+|, and pj = pT_|^+|, 7 = 
I a"*" I -I- 1 , . . . , I a+ 1 -I- 1 P~ I . The definitions are illustrated with an example in table 1 . 

Due to the conservation of the total spin the matrix elements Da^p ' 'p" (^1 , . . . , ^^) 
vanish if | a"*" | -|- | P~ | 7^ m. They are non-trivial only if | a"*" | -|- | P~ | = m. For this case 
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"2 -r, 



-1 



Figure 1 : The canonical contour C in the off-critical regime A > 1 (left panel) and in 
the critical regime | A| < 1 (right panel). For A > 1 the contour is a rectangle with sides 
at ±i| and iy, where yis slightly smaller than ^. For |A| < 1 the contour surrounds 

the real axis at a distance |y| slightly less than 



we suggest the integral representation 



n 



JtO, 



a;-i 



^.^i,yc2ra(l + o(co;)) ^Vi 



/t=aj+l 



det(-G((0;,^^)) 



ni<y</t<msh(^;fe-y sh((0;-(o^-ri) ■ 



(16) 



Equation (16) holds for any finite Trotter number A'^. The Trotter number enters the 
right hand side of (16) implicitly through the functions a((o), a(co) = l/a(co) and 
G(co, ^). Performing the Trotter limit for means to replace these functions by their 
respective Trotter limits. For brevity we show only the non-linear integral equation 
that determines the Trotter limit of a. 



lna(X) = -^- 



2/sh2(ri) 



d(£) sh(2ri)ln(l + a((o)) 



T rsh(;i)sh(X + 'n) Jc27iish(^-co + 'n)sh(A.-a)-ri) 



(17) 



The corresponding finite Trotter number equation can be found in [7]. The func- 
tion G(to, ^) which was introduced in [7] has to be calculated from the linear integral 
equation 
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sh(ri) 

sh(^-X)sh(^-X+'n) 



+ 



d(0 



sh(2ri)G(a),^) 



c 27ci(l + o(a))) sh(X-co + Ti)sh(?i-a)-ri) 



(18) 



and generalizes the 'density function' p(A,), that determines the ground state energy 
and the zero temperature magnetization, to finite temperatures and to the inhomoge- 
neous case [7]. The canonical contour C in (16)-(18) depends on r| and is shown in 
figure 1. 

Our conjecture (16) is based on the following observations, (i) Equation (16) is true 
for m = 1 , 2 which can easily be verified using the commutation relations comprised in 
the Yang-Baxter algebra and equations (98) and (104) of [7]. (ii) In two cases equation 
(16) reduces to the approximant to the emptiness formation probability, namely for 

= Py = 1 for j — 1 , . . . , m and ay = |3y = 2 for 7 = 1,..., m. These two cases 
correspond to taking the limits (p — > ±00 in our formulae for the generating function 
of the S^-S^ correlation functions in [7]. One arrives at (16) by applying the result of 
appendix C of [12] (for more details see [6]). (iii) In the zero temperature limit (16) 
reduces to the result of [14] which, in turns, generalizes the formulae [8, 9] of Jimbo 
et al. to finite values of the external magnetic field. Some important intermediate steps 
of a proof of (16) in the general case are sketched in the next section. 

Taking the Trotter limit and the homogeneous limit (for the latter compare [14]) in 
(16) we obtain an integral formula for the density matrix elements (14), 



d(£>; 



n 



,.=l^||+iJc2m(l + a(coy)) 

af-^)G(tOy,^)|^=o 



det 



sh«^-i(coy-'n)sh"'-«^(tOy) 

shP7-i((o^. + n)sh'"-P7((o^.) 

1 



{k-l)\ 



Ul<j<k<mM(^j-(^k-^) 



(19) 



We should mention that the functions 1/(1 + ci((o)) and 1/(1 + a(co)), respectively, 
appear quite naturally here, since they generalize the Fermi functions for particles and 
holes to the interacting case. 



4 Elements of a proof of the integral formula 

By analogy with the example of the generating function of the S^-S^ correlation func- 
tions treated in [7] we expect that a proof of the integral formula (16) for the gen- 
eral density matrix element can be achieved in three steps. Step 1 consists in calcu- 
lating the action of a string of operators T^^^ (^1) • • • T^"'{b,m) on the left Bethe state 
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({X,}| = (0|C(A,^/2) • • • CQ^i)- As it is clear from the structure of the Yang-Baxter alge- 
bra, the result is a linear combination of vectors of the same form with some of the "kj 
replaced with inhomogeneity parameters 1,^- Step 2 is to calculate the scalar product 
of the vectors occurring in the sum of step 1 with the Bethe state to divide by 

the norm and by the product 0^=1 Ao(^j) of transfer matrix eigenvalues, and to rewrite 
the resulting expression in a form that is suitable for taking the Trotter limit. This step 
is the same as in the derivation of the integral formula for the generating function of 
the S^-S^ correlation function in [7]. Thus, we can use our former results. 



Lemma 1. [7] 

{{^+}U{X-}\{X}) 

■ri 1 

n 



\x-\\x+\ 



sh(Xi 



^ K) 



dctG{X-f,^+). (20) 



Here we divided the sets {X} and {^} into disjoint subsets, {X^} and {^^}, and em- 
ployed the notation \X\for the number of elements in a set {X}. The function a{X) on 
the right hand side is the vacuum eigenvalue of the monodromy matrix element A{X). 

In step 3 we have to transform the sums obtained in steps 1 and 2 into integrals over 
the canonical contour C. This involves the calculation of residua and the resummation 
of the many terms emerging in this procedure. 

So far we could not complete step 3 in the general case, but only for examples of 
small m. Yet, we have obtained a complete and satisfactory result for step 1 that can 
be summarized in the following 

Lemma 2. Multiple action of monodromy matrix elements on a general state. Let 
A,i , . . . , A,M+m e C Z?^ mutually distinct. Then 

M 



■k=\ 



(o| Y[c{Xk) T^;{XM+i)...T^;;;{XM+m) = 



M+Yi M+y2 

E E 



M+Jm 

E 



Yla{Xi )c{XM^a+,Xi ) 



M+aJ 

n 

k=i 



b{Xk,Xi 



■m-\a+\ M+^j 
. j=l i j ■' k=l 



1 



b(Xe Xk) 



i- M+m 

k=l 



(21) 
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Here a{X) and d{X) are the vacuum eigenvalues ofA{X) and D{X), respectively. The 
sequences (ccj") and (p^) were defined in section 3 above. We arrange all aj and P^ 
in non-decreasing order, in such a way that Pr appears left to a+ if p^ = a+ This 
defines the sequence (y„)™^j. The position ofa^ in this sequence is denoted by yt and 
the position of^^ by (see table 1 for an example). 

Lemma 2 can be proven by induction over m. Using the fact that c(?i, X) = I, the 
well-known [18] 'elementary' commutation relations for moving A (A,), or D{X) 
through a product of Cs can be rewritten in the form 

M M+1 M+1 1 M+1 

(0| YlC{lk) A{Im+i)= <h)c{lM+iM n TrT^ {0\YlC{h), (22) 
^k=i e=i ^ t=i k=\ 

M M+1 M+1 1 M+1 

(0| ]\C(kk) D{Xm+i)^ £ J(;i£)c(;i^Xm+i) H un \ \ (oi n<^(^'t)' (^3) 

k^£ k^t 
|- M -, M+1 M+1 i-M+1 1 

(0| na^fc) B{Xm+i)^ £ £ d{Xi,)c{Xi„XM+i) Yl un 7\ 
'-k=i ^1=1 £2=1 k=i '^y^h^^kj- 



a{Xt^)c{XM+\M2) n 



M+1 



1 



k^eiA 



M+1 

(0| n Ci'^k)- (24) 

k=i 
k^iiA 



These formula prove (21) for m = 1. They can also be used in the induction step from 

n > 1 to n + 1 . 

In order to calculate the algebraic Bethe ansatz expression for the finite Trotter 
number approximantD^v^ equation(15), we have to setM = A^/2in(21), then insert the 
solution {Xjj^^j of the Bethe ansatz equation which belongs to the leading eigenvalue, 

then multiply by from the right and finally divide by {{l}\{'k})llJ=i-^o{^j)- 
The sums on the right hand side of the resulting formula are naturally divided into 
a part from 1 to M over Bethe roots and a part from M + 1 to Jj over inhomo- 
geneities The sums over the Bethe roots (and only these) can be transformed into 
integrals over the canonical contour C, since for any function /(co) holomorphic on 
and inside C the identity 



d(of{(o) 



c 27ti(l + a(co)) 



d(of{(o) 



c 27ti(l + o(to)) 



(25) 



holds. If /(co) is only meromorphic additional contributions appear due to the poles 
of /(to). In any case it is clear that all sums over Bethe roots in the expression for 
may be transformed into integrals over the contour C and that, after performing this 
transformation, a sum over multiple integrals will remain. These integrals will involve 
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not more than m integrations, and the leading term involving precisely m integrations 
will be generated by the m-fold sum from 1 to M over the Bethe roots. 



M M 

®=E E 

£l = l t2 = \ 



M 



M+aj 



.j=i y ' y 

M+P7 



1 



1. •••,«„+ 



m^J{X-}\{X}) 



k^h,...i- y 



where {X^} is the complement of in {k} = {X;}^^^, and {k'^} = {'k(,^ A,£^}. 
Inserting (20) into the right hand side of (26) and using the Bethe ansatz equations 
(see [7]) we arrive at 



M 

&= 



det(-G(^,^,)) 



£u.^m=l ni<;<fe<m M^k - ^j) sh(X4 - Xi. - Tl) 

1 



a 



k=\^'M 

m 1 

TT L_ 

.;fe=|a+|+l '^'(^4) 



|(x+|a|-l 



.j=\ k=\ 



k=a-j+l 



^-|a+|Py-l m 

n n ^h(Vi..-^^+^) n ^mv,..-^^) 

7=1 k=\ fc=P7+l 



• (27) 



Using (25) this expression turns into (16) plus a sum over terms involving less than 
m integrals. In other words, the right hand side of (16) is the unique leading term as 
described above. Thus, our conjecture (16) means that the subleading terms mutually 
cancel each other. As mentioned above we verified this statement with a number of 
examples. 



5 Discussion 

We have presented new integral (19) formulae for the density matrix of the XXZ chain 
and for its inhomogeneous generalization (16). These formulae have been verified for 
the general matrix element for small m and for the special case of the emptiness forma- 
tion probabilities for all m (note that {e\\. . . e,n\)TM ^^id {e\\. . . e^^jM are different if 
h ^ 0). We also outlined two important steps of the proof of the general formula on 
which we are working now. 

We think that our formulae have a great potential for applications in the calculation 
of finite temperature static correlation functions of the XXZ chain. We believe that the 
density matrix elements may efficiently be summed up [13]. We hope that for short 
segments the multiple integrals may be reduced to single integrals as e.g. in [10]. We 
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have started to evaluate some of the integrals numerically. Last but not least we are 
very curious if the analysis of the zero temperature inhomogeneous case as developed 
in [1] carries over to finite temperatures. 
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